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The dynamics of a large quantum spin coupled parametrically to an optical resonator is treated in 
analogy with the motion of a cantilever in cavity optomechanics. New spin optodynamic phenomena 
are predicted, such as cavity-spin bist ability, optodynamic spin-precession frequency shifts, coherent 
amplification and damping of spin, and the spin optodynamic squeezing of light. 



o 

(N 

CO 

Oh 
-)— • 

(N 
> 

CO 

in 

00 
CO 

in 

o 
o 



X 



Cavity optomechanical systems are currently being ex- 
plored with the goal of measuring and controlling me- 
chanical objects at the quantum limit, using interactions 
with light [1 . In such systems, the position of a mechan- 
ical oscillator is coupled parametrically to the frequency 
of cavity photons. A wealth of phenomena result, in- 
cluding quantum-limited measurements [2^, mechanical 
response to photon shot noise [3], cavity cooling [5, and 
ponderomotive optical squeezing [F. 

Concurrently, spins and psuedospins coupled to elec- 
tromagnetic cavities are being researched in atomic [6], 
ionic ^, and nanofabricated systems [8, 9 , with appli- 
cations including magnetometry [10 , atomic clocks [11^ . 
and quantum information processing [6l |9l [12] . In con- 
trast to mechanical objects, spin systems are more eas- 
ily disconnected from their environment and prepared in 
quantum states, including squeezed states [11. 

In this Rapid Communication, we seek to link these 
two fields by exploiting the similarities between large- 
spin systems and harmonic oscillators [13 to construct a 
cavity spin optodynamics system in analogy to cavity op- 
tomechanics. Optomechanical phenomena map directly 
to our proposed system, resulting in spin cooling and 
amplification [14^115], nonlinear spin sensitivity and spin- 
cavity bistability [I6l [17] , and spin opto-dynamic squeez- 
ing of light [14l[T8]. Such a system may find application 
as a quantum-limited spin amplifier or as a latching spin 
detector. We detail these phenomena using currently ac- 
cessible parameters, and we propose realizations either 
using cold atoms and visible light or using cryogenic solid 
state systems and microwaves. 

An ideal cavity optomechanics system, consisting of 
a harmonic oscillator coupled linearly to a single-mode 
cavity field, obeys the Hamiltonian 



H = hujcfi -\- hWzO/a — fziio (a^ -\- a) n -\- H 



in/ out' 



(1) 



Here a is the oscillator's phonon annihilation operator, 
h is the photon number operator, uj^ is the natural fre- 
quency of the oscillator in the dark, and ujc is the bare 
cavity resonance frequency. / is the radiation-pressure 
force applied by a single photon, while zho = \/h/2muJz 
is the harmonic oscillator length for oscillator mass m. 
Hin/out describes the coupling of the cavity field to exter- 
nal light modes. Under this Hamiltonian, the cantilever 



position z and momentum p evolve as dz/dt = p/m and 
dp/dt = —muolz + ffi. 

To construct a spin analogue of this system, we con- 
sider a Fabry-Perot cavity with its axis along k (Fig. 
[l]). For the collective spin, we first consider a gas of 
N hydrogenlike atoms in a single hyperfine manifold of 
their electronic ground state, each with dimensionless 
spin s and gyromagnetic ratio 7. The atoms are op- 
tically confined at an antinode of the cavity field. An 
external magnetic field B = Bh is applied to the atoms. 
The detuning Aca between the cavity resonance is cho- 
sen to be large compared to both the natural linewidth 
and the hyperfine splitting of the atoms' excited state. 
In this limit, spontaneous emission may be ignored and 
the single-atom cavity-field interaction energy, Hstark = 
(%o/^ca)^ (1 ± ^;k • s), comprises the scalar and vector 
ac Stark shifts [19 , where qq quantifies the atom-cavity 
coupling and v the vector shift. 

Summing over all atoms q^ we obtain the system 
Hamiltonian: 

U = hujc (n+ + n_) + Uin/out + X^ f - ^tB • s^ 

q ^ 

+ ^[(n++n_) + t;(n+-n_)k-s,]Y (2) 

with number operators h± for the a^ polarized optical 
modes. 

The above Hamiltonian can be rewritten as the inter- 
action of the collective spin operator S = ^ Sg with an 
effective total magnetic field Beff = r2eff/7, giving [20] 



r^eff = l^Lb + Vtc (^+ — n_) k. 



(3) 



Here VLl = jB and ftc = —'^oil^ca- Altogether, the 
cavity spin optodynamical Hamiltonian is 



H = n CJ, 



^9l 



{n^^n-)^nin/out-^^ei^-S. (4) 



Now consider the external magnetic field to be static 
and oriented along i, orthogonal to the cavity axis. In 
the limit (S) o:^ 6'i, the spin dynamics become 

^ = nLSk-n,s{n+-h.), ^ = -nLS,. (5) 



r^c(n^-n_)k 




FIG. 1. (Color) An ensemble of atoms trapped within a driven 
optical resonator experiences an externally imposed magnetic 
field along i and a light- induced effective magnetic field along 
the cavity axis k. The evolution of the collective spin S re- 
sembles that of a cantilever in cavity optomechanics. 



The analogy between cavity optomechanics and 
spin optodynamics is established by assigning 
^ -^ -zuoSk/ASsQL and p -^ puqSj/ASsql, 
where zuo and pno = ^/(2^ho) are defined with 
ujz -^ ^L [13] and A^sql = \/5'/2 is the standard 
quantum limit for transverse spin fluctuations. Eqs. ([5| 
now match the optomechanical equations of motion 
with the optomechanical coupling defined through 
/zho^^— ^^c^'S'sql(^+— ^-)- The main result of this 
work, that various cavity optomechanical phenomena 
are manifest also in cavity spin optodynamical systems, 
is immediately established. 

Let us now elaborate on these phenomena. To obtain 
general results, we will proceed without assuming S c^ 51, 
except in certain cases, noted in the text, where some 
physical insight is gained. We begin with effects for which 
both the light field and the ensemble spin may be treated 
classically, i.e. by letting S = (S) and n± = {n±). 

Cavity-spin bistability: We start with the static 
behavior of the system by finding the fixed points of the 
system. The collective spin vector is static when S is 
parallel to r^eff. Writing S = 5'(isin^o + kcos6>o), this 
condition requires n.^ — fi- = (Ql/^c) cot^o- The intra- 
cavity photon numbers are determined also by the stan- 
dard expression for a driven cavity of half line- width /^, 
i.e. n± = nmax,±[l + (Ap,± ± l]c^cos6>o) /i^'^]~^ with 
uj± = {uJc + Ng^/Aca) + ^p,± being the frequency of 
laser light of polarization a^ driving the cavity and 
^max,± characterizing its power. These two expressions 
for n+ — n_ may admit several solutions (Fig. ^. 

As typical in instances of cavity bistability [22 , sev- 
eral of the static solutions for the intracavity intensities 
may be unstable. To identify such instabilities, we con- 
sider the torque on the collective spin when it is dis- 
placed slightly toward +k from its static orientation. 
Stable dynamics result when such displacement yields a 
torque N • j with the sign a = sgn(sin^o)- Geometri- 
cally, this stability requires that the spin vector be dis- 
placed further in the +k direction than the vector aHeff. 
Quantifying the linear response of the intracavity effec- 




■> ■ v^^ 



15 




-10 -5 5 10 -10 -5 5 

Ap/K Ap/tz 



FIG. 2. (Color) Cavity-spin bistability in a cavity driven 
with linearly polarized light. We consider N = 5000 spin- 
2 ^'^Rb atoms, Qc/k^ = 1.25 xH)"^ Ql/k = 3.3 x 10"^ and 
nmax,± = 15 (similar to Ref. 21 ). (a) As Ap is varied, several 
stable (black) and unstable (gray) static spin configurations 
are found. Configurations for A^/ai: = — 4.8 are depicted, (b) 
The cavity exhibits hysteresis as the probe is swept with posi- 
tive (dashed blue) or negative (red) frequency chirps, with the 
spin initially along i. Rapid transitions as Ap/ n is swept up- 
ward from -2.8 or downward from involve symmetry break- 
ing as the cavity becomes birefringent; we display n+ and fi- 
assuming the stable branch closer to ^o = is selected. Here, 
Aca/27r = 20 GHz from the D2 transition, ^o/27r = 15 MHz, 
A>:/27r = 1.5 MHz. 



tive magnetic field to variations of the collective spin via 
A = Vtcd{nj^ — n-)/dSk^ the static spin orientations are 
found to be unstable when a\ > Ql\ csc^ Oo\/S. 

Opto-dynamical Larmor frequency shift: The 
dynamics of the spin precessing about one of the stable 
configurations can be parameterized by the precession 
frequency, which is shifted from Ql by two effects. First, 
there is an upward frequency shift from the static modi- 
fication of the effective magnetic field, leading to preces- 
sion at the frequency Ql = ^l\ csc^oI when A = 0. A 
second shift occurs when the spin dynamics are slow com- 
pared to the response time of the cavity field {Ql ^ ^)- 
Here, the precessing spin modulates the cavity field, 
which, in turn, acts back upon the spin to modify its 
precession frequency. When the precession amplitude is 
small, a solution of the spin equations of motion derived 
from the Hamiltonian in Eq. Q yields an overall preces- 
sion frequency CIl^ where 



Qt 



Qj 



XftLSsinOo. 



(6) 



The quantity ks = —XQlS sin Oq serves as the analogue 
of the optical spring constant [23] , and leads to shifts of 
the Larmor precession frequency with a sign and mag- 
nitude that depend on the spin orientation, A and the 
frequency, intensity, and polarization of the cavity probe 
fields. When the precession amplitude is large, the dy- 
namics become essentially nonlinear. In this case the dy- 
namics can be described by numerical simulation (Fig.lsl). 
Coherent amplification and damping of spin: 



Now we consider the effects of the finite cavity response 
time K:~^ on the spin dynamics. To develop an intu- 
itive picture, we consider the unresolved sideband regime 
Ql < /^, in a frame (indicated by the index "r") corotat- 
ing with the collective spin, with ir aligned to the fixed 
point. We assume the spin to be precessing at a near 
constant rate, and the cavity field response to this pre- 
cession to be simply delayed by hz~^. Employing the 
rotating-wave approximation, the delay causes the ef- 
fective field r2eff,r to point out of the ir-k^ plane, with 
^eff,r • jr = — (ttAS'/cr siu^ ^q siu 0)/2, whcrc (j) = Ql'/k,. 
The collective spin now experiences a torque in the kr 
direction, giving 



(a) 5000 



dSk,r —aX sin^ Oq sin (l)Si^r 
dt ~ 2 



S^ 



k,r 



(7) 



For positive (negative) values of aX , the Larmor pre- 
cession frequency is shifted down (up) and the spin is 
damped toward (amplified away from) its stable point. 
Similar relations apply to cavity optomechanics ^ . The 
deflection of the spin toward or away from the stable 
points persists for large precession amplitudes (Fig. |3|. 

This cavity-induced spin amplification or damping dif- 
fers from conventional optical pumping in two important 
respects. First, while the spin polarization generated by 
optical pumping relies on the polarization of the pump 
light, the target state for cavity-induced spin damping 
is selected energetically. Similar to cavity optomechani- 
cal cooling [4j, cavity enhancement of Raman scattered 
light drives spins to the high- or low-energy spin state 
according to the detuning of probe light from the cavity 
resonance, independent of the polarization. Second, this 
amplification or damping of the intracavity spin is coher- 
ent, preserving the phase of Larmor precession, at least 
within the limits of a quantum amplifier. 

Spin optodynamical squeezing of light: We now 
consider quantum optical effects of cavity spin optody- 
namics. One such effect is the disturbance of the collec- 
tive spin due to quantum optical fluctuations of the cav- 
ity fields. In cavity optomechanics, intracavity photon 
number fluctuations disturb the motion of a cantilever, 
providing the necessary backaction of a quantum mea- 
surement of position [25] . The analogous disturbance of 
optically probed atomic spins (or pseudo-spins) has been 
studied both in free-space [26] and intracavity [11] [27] 
configurations. In an optomechanics-like configuration, 
e.g. with B (X i, backaction heating of the atomic spin 
enforces quantum limits to measurement of the precess- 
ing ensemble and also set limits on optodynamical cool- 
ing. In contrast with optomechanical systems, optically 
probed spin ensembles readily present the opportunity to 
perform quantum- non-demolition (QND) measurements; 
with B oc k, the detected spin component Sk is a QND 
variable representing the energy of the spin system. 

The noise-perturbed spin acts back upon the cavity 
optical field, mediating a self-interaction of the light field 
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FIG. 3. (Color) Simulations of spin dynamics for S = 5000, 
Ql/27v = 200 kHz, Qc/27r = -2.3 kHz, At/27r = 1.8 MHz, 
n+ = 10, and Ap,+ = 0.37ai:. (a) Time evolution of Si (black) 
and Sk (blue), following spin preparation near i, shows ampli- 
fication, reorientation, and damping toward the high-energy 
stable orientation near — i. Note the different scales on the 
horizontal axis, (b) Logarithmic optical spectral noise power 
relative to that of coherent light, plotted vs. quadrature an- 
gle (/) (amplitude quadrature at = 0), shows inhomogeneous 
optical squeezing. Simulation results shown in color, and lin- 
earized theory (Eq. 10) as contour lines every 5 dB. 



that can result in optical squeezing. To exhibit this ef- 
fect, we consider a cavity illumined with (j+ circular po- 
larized probe light with detuning Ap. The dynamics of 
the cavity field are given by 



dc-^ 



= {iAp - /^ + incSk)c^ -^i^(r]^ ^V) . (8) 



Here, r] gives the coherent-state amplitude of the drive 
field and the noise operator ^+ represents its fluctua- 
tions. When evaluating the dynamics numerically, we 
consider a semiclassical Langevin equation, converting 
^+ into a Gaussian stochastic variable with statistics re- 
lated to those of the noise operator, and replacing the 
operators c+ and S with c-numbers. This substitution is 
appropriate for moderately large values of n and S. 

Fig. [3] portrays the simulated evolution of a spin pre- 
pared initially in a low-energy spin orientation (close 
to i), driven by a blue-detuned cavity probe. Coher- 
ent spin amplification directs the spin toward the stable 
high-energy configuration (near — i), yielding a dynami- 
cal steady state characterized by a negative temperature. 

To obtain analytical expressions for the evolution dy- 
namics, we follow the example of cavity optomechanics |5| 



by linearizing the Langevin equations for spin and opti- 
cal fluctuations about their steady-state value. The spin 
projection Sk responds to amplitude-quadrature fluctua- 
tions of the cavity field £,a{^) with susceptibility 



xH 



—Ql ^c^/^-\- 



UH Ql'^ + ksR{uj) -uj'^^ iujToiuj) 



(9) 



^i 



where To{uj) = 2/^ ^2+^2 
spin damping, and R{uj) 



is the cavity optodynamic 
. The susceptibil- 



ity is largest for uj c^ O^l" - The driven spin feeds the 
fluctuations back onto the cavity field, yielding the in- 
tracavity field fluctuation spectrum 



Qr 



ksR{uj) 



■ iujTo{oj) 



(10) 

where £,p{cjo) is the input spectrum of phase fluctuations. 
This fluctuation spectrum exhibits inhomogeneous opti- 
cal squeezing (Fig. [SJd). 

Applications: The analogy of cavity optodynamics 
widens the range of phenomena accessed through the ma- 
nipulation and detection of quantum spins within opti- 
cal cavities, enabling several applications. For example, 
bistability in cavity-coupled single-spin systems serves to 
increase the readout fidelity of cavity-coupled qubits [28] . 
Similarly here, cavity-spin bistability could be used as a 
Schmitt trigger for the collective spin: if the probe power 
is turned on diabatically in the bistable regime, the cav- 
ity transmission will latch into either a bright or dark 
state, depending on whether the initial spin state is be- 
low or above a parametrically chosen threshold value. 
The cavity-spin system may also be used as a phase- 
preserving amplifier for spin dynamics occurring near 
the shifted precession frequency ^l\ with amplification 



noise given in Eq. (10). Both applications may aid mea- 
surements of ac magnetic fields, amplifying weak signals 
above technical sensitivity limits. 

Conversely, cavity spin optodynamics may be applied 
as a powerful simulator of cavity optomechanics, with the 
spin system allowing for new means of control. For exam- 
ple, precession frequencies may be tuned rapidly by vary- 
ing the applied magnetic field, simulating optomechanics 
with a dynamically variable mechanical spring constant. 
Alternately, spatial control of inhomogeneous magnetic 
fields may be used to divide a spin ensemble into sev- 
eral independent subensembles, simulating optomechan- 
ics with several mechanical modes. 

In addition to the dilute gas implementation discussed 
so far, a similar system could be constructed using solid- 
state spin ensembles and microwave resonators. For ex- 
ample, using ensembles of nitrogen- vacancy defects in di- 
amond coupled to the circular polarized evanescent radi- 
ation of a crossed microwave resonator [29 , the Hamilto- 
nian of Eq. Q is obtained by the ground rris = ±1 elec- 



tronic states as the pseudospin and replacing the ac Stark 
shift with an ac Zeeman shift from microwave radiation 
near the 2. 8- GHz crystal- field-split Zeeman transition. 
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